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PACS 05. 40. -a - Fluctuation phenomena, random processes, noise, and Brownian motion 
PACS 02. 50. -r - Probability theory, stochastic processes, and statistics 
PACS 75. 10. Hr - Spin-glass and other random models 

Abstract. - We study the joint probability distribution function (pdf) Pt(M, Xm) of the maxi- 
mum M of the height and its position Xm of a curved growing interface belonging to the univer- 
sality class described by the Kardar-Parisi-Zhang equation in 1 + 1 dimensions, in the long time 
t limit. We obtain exact results for the related problem of p non-intersecting Brownian bridges 
where we compute the joint pdf P P {M, tm), for any finite p, where tm is the time at which the 
maximal height M is reached. This yields an approximation of Pt(M, Xm) for the interface prob- 
lem, whose accuracy is systematically improved as p is increased, becoming exact for p — > oo. We 
show that our results, for moderate values of p ~ 10, describe accurately our numerical data of a 
prototype of these systems, the polynuclear growth model in droplet geometry. We also discuss 
applications of our results to the ground state configuration of the directed polymer in a random 
medium with one fixed endpoint. 
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Introduction. — The study of fluctuations in stochas- 
tic growth processes has attracted much attention during 
the last two decades [TJ [2] . Such processes are ubiqui- 
tous in nature as they appear in various physical situ- 
ations ranging from paper wetting to burning fronts or 
growing bacterial colonies. In many experimental settings 
the growth starts from a point like region (a seed) with a 
strong tendency to evolve towards the approximate spher- 
ical symmetry. Such examples include fluid flow in porous 
media, adatoms and vacancy islands on surfaces [3| but 
also biological systems such as tumors 

To describe such phenomena driven by a growing in- 
terface, several models have been studied, like the Eden 
model, polynuclear growth models (PNG) or ballistic de- 
position models, among others [2j. In 1 + 1 dimensions, 
it is widely believed that all these models belong to the 
same universality class as that of the Kardar-Parisi-Zhang 
(KPZ) equation [5|[6]. At time t, the width of the in- 
terface Wl(£), for such a system of size L, behaves like 
Wi,{t) ~ L^W^/ L z ) with universal exponents C — \ and 
z = | 7|. In the growth regime to <C t <C L z (where to is 
a microscopic time scale) exact results for different lattice 
models have shown that the notion of universality extends 
far beyond the exponents and also applies to full distribu- 
tion functions of physical observables [8]- 10 . In particular, 



with ft = 2 (respectively j3 = 1) for the curved geometry 
(respectively for the flat one), which describes the edge 
of the spectrum of random matrices in the Gaussian Uni- 
tary Ensemble (respectively of the Gaussian Orthogonal 
Ensemble) [ll]. Height fluctuations were measured in ex- 
periments, both in planar |12| and more recently in curved 
geometry in the electroconvection of nematic liquid crys- 
tals [13] and a very good agreement with TW distributions 
was found. 

Here we focus on a prototype of these models, the PNG 
model |14| , but our results hold more generally for curved 
growing interfaces in the KPZ universality class (see be- 
low). It is defined as follows. At time t — a single island 
starts spreading on a flat substrate at the origin x = 
with unit velocity. Seeds of negligible size then nucleate 
randomly at a constant rate p — 2 per unit length and 
unit time and then grow laterally also at unit velocity. 
When two islands on the same layer meet they coalesce. 
Meanwhile, nucleation continuously generates additional 
layers and in the droplet geometry, nucleations only oc- 
cur above previously formed layers. Denoting by h(x,t) 
the height of the interface at point x and time t, one thus 
has h(x,t) = for |a;| > t. On the other hand, in the 
long time limit, the profile for \x\ < t becomes circular 
h(x, t) 



2ty/l - (x/ty 



the scaled cumulative distribution of the height field co- 
incides with the Tracy- Widom (TW) distribution 7^(0 
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but there remain height fluc- 
tuations around this semi-circle. A natural way to charac- 
terize these fluctuations is to consider the maximal height 
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Fig. 1: a): Height profile h(x,t) for fixed time t for the PNG 
model, b): An optimal path of the DPRM of length t. c): A 
p- watermelon configuration (p = 8). d): Sketch of the method 
to compute the joint pdf P f (M, tm). 



M and its position Xm (see Fig. 
KPZ scaling, one expects M — It 



[I| a)). According to 
~ tz, while Xm ~ £ s 
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The purpose of this Letter is to provide an analytic 
approach to the joint pdf Pt(M, X M ) for large t. 

The PNG model can be mapped onto the directed poly- 
mer in a random medium (DPRM) on the square lattice 
with one fixed end 2fl0 . In this language M is the ground 
state energy while Xm is the transverse coordinate of the 
free end of the optimal polymer of length t (see Fig. [I] 
b)). Related questions for the continuum DPRM are cur- 
rently under active investigations 17 . The marginal dis- 



tributions of M and Xm are already interesting and re- 
lated extreme value quantities, like the maximal relative 
height, have been extensively studied in the stationary 
regime t ^> L z 18-20 . Much less is known in the growth 
regime which we focus on. From the mapping onto the 
DPRM, one identifies the pdf of the maximal height M in 
the droplet geometry with the pdf of the height hn a t(x,t) 
at a given point x and time t in the planar geometry [21] . 
Therefore we conclude that the pdf of M, suitably rescaled 
and shifted, is given by [l6j[22j, the TW distribu- 

tion for = 1. On the other hand, the computation of 
the (marginal) distribution of Xm is a challenging open 
problem |16| . 

To compute P t (M, Xm), we exploit the exact mapping 
between the height field in the PNG model and the top 
path of p non-intersecting random walkers, so called vi- 
cious walkers 23 1 , in the limit p — > oo [8 24 . Here we 
consider "watermelons" (Fig. [l]c)) where p non-colliding 
Brownian motions X\(t) < ■ ■ ■ < x p (r) on the unit time 
interval are constrained to start and end at (i.e. Brow- 
nian bridges). In the large p limit, one can show, using 
the connection between this vicious walkers problem and 
random matrix theory, that x p {r) also reaches a circular 
shape of amplitude y/p, x p (t) ~ 2^/py/ r(l — r) while the 



fluctuations are in that case of order p~® [25]. Hence x 
and t map onto h and x in the growth model while p plays 
the role of t. This mapping, for p,t > 1 reads [8 24 



h{utl,t)-2t _ x v {\ + l jp 2 

1 = n = A2W — u {i) 

13 p e 

where A2(u) is the Airy2 process [8] which is a sta- 
tionary, and non-Markovian, process. In particular, 
ProbaL42(0) < £] = -7"2(£)- ^ n tms Letter, we compute 
exactly the joint distribution P p (M,tm) of the maximal 
height M and its position tm for the vicious walker prob- 
lem (Fig. [l]c)). While the maximal height M has been 
recently studied 25 26 , nothing is known about the dis- 



tribution of tm , which has by the way recently attracted 
much interest in various other one-dimensional stochastic 
processes [27}]30] . Our results are not only relevant, for fi- 
nite p, for the vicious walkers problem, but thanks to the 
above relation ([!]), become exact for p — > oo, for curved 
growing interface, as well as for the DPRM. We actually 
show that for moderate values of p ~ 10, our analytical 
formula describes quite accurately our numerical data for 
the PNG model. 

Method. — The basic idea of the method to compute 
P P (M, tm) is to divide the "watermelons" configuration in 
two time intervals, r £ [0, tm] and r £ [tm, 1] and use the 
Markov property of the whole process to treat these two in- 
tervals independently (Fig. [l]d)). In both intervals, the p 
vicious walkers are constrained to stay below M, while we 
impose x p (tm) = M. To compute the propagator of these 
constrained vicious walkers in each sub-interval, we use 
a path-integral approach. Let us denote p<jvf(b, ib|a, t a ) 
the propagator of p non-intersecting Brownian motions, 
starting in a e (a 1; -- - ,a p ) at time t a and ending in 
b = (&].,•■• , b p ) at time tb and constrained to stay below 
M in the time interval [t a ,tb\. p<jvf(b, i&|a, t a ) is given 

by the sum of the weights exp — \Y?h=i It* (^r) 2< ^ T 
over all trajectories satisfying X\(t) < x 2 (r) < ••• < 
x p (t) < M, for t £ [t a ,tb\. In the language of path in- 
tegrals, this correponds to the propagator (in imaginary 
time) of p quantum free fermions with an infinite wall 
in x = M, the associated Schrddinger Hamiltonian being 
Hm = — 2 ELi ®Xi + Vm(x). The hard wall potential is 
given by V M (x) = if x < M and V M (x) = +oo if x > M 
|25| . The use of fermions incorporates naturally the non- 
colliding condition 25 31 32 . This allows to write this 
propagator as 



P<M (b,t b \ a ,t a ) = (b|e-( t! >-*^|a) 



(2) 



which we can compute using a spectral decomposition over 
the fermionic eigenfunctions of Hm- Before doing this, 
we notice that the "watermelons" configurations that we 
study here are actually ill-defined for Brownian motions 
which are continuous both in space and time. It is in- 
deed well known that if two walkers cross each other they 
will recross each inifinetly many times immediately after 
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the first crossing. This means in particular that it is im- 
possible to impose Xj(0) = Xj+i(0) and simultaneously 
x,(0+) < x i+ i(0+). Here, following Ref. [l9j[25j[27] , we 
adopt a regularization scheme where we impose that the 
p walkers start and end at < e < •• • < (p— l)e and take 
eventually the limit e — > 0. We use an additional cut-off 
procedure by imposing that x p (tm) = M — r\ and then 
take the limit 77 — 5- 0. 

Results for p vicious walkers. — The calculation of 
P p {M,tm) requires the computation of p<Af(b, t(,|a, t a ). 
Expanding Eq. ([2| over the fermionic eigenvectors of Hm 
yields 

p <M (b,t b \a,t a ) = (W (tb ~ ta)HM \a) 



dk e-V (*»-*.) det fatbi) det c/>* k ( aj ) , (3) 



where <j>h{x) — y f sin [k(M — x)\ naturally appear as 
the eigenvectors of Hm and where we use the notations 

S' J dk = f °° dki ■ ■ ■ / °° dk p and k 2 = k\ + ■ ■ ■ + k 2 p . In Eq. 
, the determinants appear as Slater determinants in the 
associated fermions problem. From this propagator ^ we 
compute P p (M,tm) by dividing the configuration in two 
time independent intervals, r £ [0,tm] and r £ [tm, 1] as 
explained above (see also Fig. [l]d)) to obtain: 

M-n 

P p (M,t m )= inn — / dyp< M (e, l|y,r Af ) 

—00 

xp<M(y,T M \e,0)S(y p - (M-tj)) , (4) 

where the delta function enforces x p (tm) = M — T) and 
where the amplitude Z p , which depends explicitly on 
e and rj is determined by the normalization condition 
Jq 1 " 00 AM J q 6.tmP p {M,tm) — 1< Using the above formula 
for the constrained propagator ^ in Eq. Q, and taking 
the limits e, rj — ¥ one obtains 



P p (M,t m ) = z- 1 M-^ 2+ 



3) 



dq / dq'e 



2M 2 



x e 



q P q' p @ P (qi,--- ,q P )® P (qi,--- ,%), (5) 



da; det fi (xj 



det g t ( Xj ) 



! det 

1<2,J<P 



dxfi(x)gj(x), (7) 



for any suitable functions /i(x) and gi(x). These integrals 
can be expressed in terms of Hermite polynomials H n (x) 
and we obtain finally 

P p {M, t m ) = B p [det D] t U(r M )D- 1 U(l - t m ) , (8) 

with B- 1 = (2tt)5 H p .zl (j!2 J ) and where D = D(M) is a 
p x p matrix 

Dij = (-l) 4 - 1 i/ 4+J - 2 (0) - e- 2M2 i7 l+ ,_ 2 (y2Af) , (9) 
while U(r) = U(M, r) is a column vector given by 

Ui(r) - r"^ fT< (m/V2t) e"^. (10) 




Fig. 2: Contour plot of P P (M, r M ) for p - 
correspond to 0.1, 0.4, • • ■ , 2.2. 



3. The contour lines 



where O p (q) is the following determinant 

6 p (q) = det q{ _1 cos - jn/2) . (6) » = 

i<»>i<p 

To compute z p , we use that P p (M, tm)Atm must yield 
back the expression for the distribution of the maxi- 
mum as computed in Ref. [25]. This allows to obtain 

z p = 7r 1+ 2 2~^ Ilj=o J- ^ n -^l- |^}> one can ex pand the 
determinant by minors and then perform the integrals over 
qi using the Cauchy-Binet identity which reads 



In Fig. Q, we show a contour plot of P p (M,tm) f° r 
p = 3. For fixed tm € [0, 1], P p (M, tm), as a function of 
M has a simple bell shape. Its behavior for fixed value of 
M as a function of tm is more interesting. For sufficiently 
large M > M*, it has a bell shape, with a maximum in 
t m = § , while for M < M* , it has a " M-shape" with two 
distinct maxima, tm — \ being a local minimum. One 
observes that M* is a slowly increasing function of p. 

By integrating our expression in Eq. (JsJ over tm, °ne 
checks that we recover the formula for the pdf of M, as ob- 
tained in Ref. [26|. Indeed, one finds Proba[max r x p (r) < 
M] = detD/nf=o (j' !2J ), where the matrix D is defined 
in Eq. ([9|. On the other hand, by integrating Eq. ^ 
over M, one obtains an expression for P p (tm)- While for 
p = 1, P\{tm) = 1, one obtains, for instance, for p = 2 



%) = 4 1- 



1 + 10t m (1-tm) 
(1 + 4t m (1-tm)) 5 / 2 



(11) 



For generic P p {tm) is a function of tm(1 — tjw) with the 
asymptotic behavior, for small tm 

P p {tm) - T^ p \ with u(p) = (p 2 +p- 2)j2. (12) 

Note that for p independent Brownian motions (without 
the non-crossing condition) one has -Ppjrce^ju) ~ t m~ 1 
29 so that the exponent u(p) bears the signature of the 
non-colliding condition. 
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The above approach can be extended to study the ex- 
treme statistics of p non-intersecting excursions, i.e. vi- 
cious walkers starting and terminating at the origin but 
with the additional constraint that they all stay positive 
inbetween |25[|33[[34| . The slight modification in our com- 
putation is to replace the Hamiltonian Hm (with a wall in 
x = M) with the box Hamiltonian Hb ox (with two walls: 
in x = and in x — M). Hence the energy levels are 
discrete, and consequently one has to treat discrete sums 
instead of integrals as before. Using the appropriate prop- 
agator in Eq. Q, one finds, taking the limits e, r\ — > 0, the 
joint pdf for the non-intersecting excursions 



E 

ni_,...,n p ,n' >0 
p-1 

x ~[[ n i n l n p \( n l, ■ ■ ■ , n l)\( n l, ■ ■ ■ ,n' p ), (13) 
t=l 

where we use the notation A p (Ai, . . . , X p ) = 
rii<i<i<j)(A? — ^») f° r tne Vandermonde deter- 
minant , and with the normalization z Pi e = 

2 S rT(p)U P j Zl(2 : > 'j!r(| + ]))Tr-( 2p2+p+2 \ Performing 
similar manipulations as before, one finds the joint pdf 
of the couple (M, tm) for the excursions configuration 
in a dctcrminantal form, reminiscent of the watermelons 
case Q: 

PpM^^m) = C P [detD B ] *U £ (r M )D £ 1 U B (l - r M ) , 

(14) 

with C- 1 = i-iy+^P 2 ' 1 ^ JYjZlU^il + The 
p x p matrix = T)e(M) is (for 1 < i,j < p) 



D ^ = E % +3 -i)(V2M n )e- 2M ", (15) 

n— — oo 

and appears in the expression of the cumulative distri- 
bution of the maximum |33| , Proba[max r x p (r) < M] = 
(-l) p detD B /(2P 2 Yl? =1 (2j-l)l). The vector elements are 
(for i = l,...,p) 



Application to stochastic growth processes. 

We now come back to the joint pdf Pt(M, Xm) for curved 
growing interfaces (Fig. [T]a)). From Eq. ^ one obtains 

P t (M,X M ) ~ i- x P Airy ((M - 2t)t-$,X M r$) (18) 

where V^hy (y, %) is the joint distribution of the maximum 
y and its position x for the process A2(u) — it 2 , u £ K. 
From Ref. [16 2Tj 22 one obtains that the (marginal) pdf 
of y, VAiiy,M(y), is given by TW, VMryMiv) = ^i(y)> 
and hence P t (M) ~ t~aJ^((M -2t)/ti). On the other 
hand, from Eq. ([I]), T'Ai T y{y,x) can be obtained from our 
results for P p (M, tm) in Eq. in the scaling limit 

y = (M — ^/p)p^ as well as x = 2{tm — \)p^ fixed one 
has indeed P p (M,tm) ~ 2p^Vxity{y,x). While our ex- 
pression ^ should be amenable to an asymptotic analy- 
sis for large p, yielding an explicit exact expression for 
PAiTyiVi x )i such an analysis deserves further investiga- 
tions. Here we show instead that, for finite values of p, 
this expression ^ describes quite accurately our numer- 
ical data for the extreme statistics of the PNG model in 
the droplet geometry. To illustrate this, we have computed 
numerically the (marginal) distribution of the position of 
the maximum Xm- We recall that, for the DPRM as in 
Fig. ([I] b)), Xm is the position of the free end of the 
optimal polymer. 
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(line), 2p 1/3 (x„ 



1/2) (dots) 
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M 
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' M 2 



71=1 



(~1) 



1 ) n e 



(16) 



The integration of formula ( 13 1 with respect to M gives 



the pdf of the time to reach the maximum in this excursion 
process. For p = 1, the formula in Eq. (13 1 yields back 
the result obtained in Ref. 127]. For p — 2, we give the 
expression of P2,e(t~m) 



Pi,e{tm) = a 



(-ir^n\nlnl{n\-nl){nl-nl) 



ni>0 



(n? 



TMnl + (1 - r M )-T-i) e 



with a = 1280/vr. 



(17) 



Fig. 3: On linear scales, the solid line is a plot of t3 P t (XM) as 
a function of Xu/tz , for t = 768, computed numerically for 
the PNG, while the dots correspond to our analytical results 
P p (r Af )(2p3) Q as a function of2ps(r M -§) forp= 10,12,14 
and p = 16. 

In Fig. [3j we show a plot of the rescaled distribution 
t3 P t [X.M A s ) as a function of the rescaled variable Xm/* 3 
for t = 768 in solid line. We also plot our exact analytical 
results for watermelons, i.e. P p {tm) / {2p~s) as a function 
of (2p3(rjvf — |)) for p — 10, 12, 14 and 16 which is com- 
puted from Eq. S as P p {t m ) = J °° dMP p (M,T M ). We 
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emphasize that the good collapse of the different curves is 
obtained without any fitting parameter. 
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Fig. 4: Log-linear plot of: the rescaled distribution of the po- 
sition of the maximum for PNG (solid line), the rescaled dis- 
tribution of the time to reach the maximum for the p = 16 wa- 
termelons configuration (dots) , and the closest Gaussian curve 
(dashed line), that do not match the PNG curve. The scaled 
quantities appears in the x, y labels (and are detailed in the 
caption of Fig [3]l . 

In Fig. |4j we show a plot of the same quantities as in 
Fig. [3] in log-linear plot. As expected, one observes some 
deviations between the result for the PNG and our com- 
putation for finite p = 16 in the tail of the distribution. 
Besides, our numerical data suggest that the marginal dis- 
tribution of Xm is non- Gaussian (the best Gaussian ap- 
proximation being shown as a dotted line). Instead, for 
Xm 3> t 2 / 3 , our data are compatible with a stretched ex- 

ponential behavior ~ t se 1 \~t s TS ) with 5 ~ 2.5 

while its precise determination requires more numerical 
efforts. Interestingly, for p free Brownian bridges, one 
can show [29] that the distribution of t m converges to a 
Gaussian distribution centered in \ of width (81og7i)~5 . 
Therefore, the non-Gaussianity of P t {X M ) is a clear sig- 
nature of the correlations in the associated vicious walkers 
problem. 

The results obtained in the present Letter extend far 
beyond the PNG model. Indeed our results hold more 
generally for stochastic growth models, and physical sit- 
uations, where the interface h(x,t) at point x and time 
t evolves according to the one-dimensional KPZ equation 



h(x,t) = uV 2 h{x,t) + ^(Vh(x,t)) 2 + C(x,t) , (19) 



di 

in a curved geometry. In Eq. ( [l9| , C(x,t) is a 
Gaussian white noise of zero mean and correlations 
(((x,t)((x',t')) = DS(x - x')S(t - t'). In Ref. (8j, it 
was shown that the fluctuations of physical obervables 



of the height field, evolving according to Eq. (19), are 
universal, up to two non-universal parameters A and 
A = D/2v. Once A and A are fixed, these fluctu- 
ations are characterized by universal distribution func- 
tions. In the curved geometry, A can be simply mea- 
sured as the radial growth rate, (h(x — 0, t)) ~ Xt 
while A can be extracted from the width of the interface 
((h(0,t) - (/i(0,i))) 2 ) 1/2 - AL/6 for a system of finite 
size L [2l]. For instance, the fluctuations of the height 
field are given by h(0,t) ~ A< + (A 2 Ai/2) 1/3 x G UE where 
Xgue is a random variable distributed according to the 
TW distribution T% (8j. This universality was convinc- 



ingly demonstrated in recent experiments 13 . There- 
fore, for an interface described by the KPZ equation ( 19 1 
in a droplet geometry, one expects that the distribu- 
tion of the position of the maximum Xm will have the 
scaling form P t {X M ) ~ CW'^Airy.Xw {X M /£(t)) with 
£(t) = (2/A)(A 2 \t/2) 2 / 3 is the correlation length. The 
scaling function VAiry,x M ( x ) i s universal and can be com- 
puted from our formulas in Eq. (|8 18 1, see also Fig. [3] 



Conclusion. — To conclude, we have obtained an ex- 
act expression for the joint distribution of the maximum 
M and the time tm at which this maximum is reached for p 
non-intersecting Brownian bridges ([8| and excursions ( 14 ) . 
We have shown that our analytic expression for moderate 
values of p Brownian bridges describe very accurately the 
extreme statistics of curved growing interfaces, becoming 
eventually exact in the limit p — > 00. In addition to their 
relevance to the DPRM, our results may have applications 
to various other situations like step fluctuations in faceted 



crystals [35] or dimer-covering problems 36 , where it was 
shown that the fluctuations are governed by the very same 
process .42 (it) — u 2 that we have studied here. Finally, 



in view of recent progresses 12 13 , it seems possible to 



observe these extremal statistics for curved growing inter- 
faces in experimental situations like nematic liquid crystals 
131. 



We thank T. Sasamoto for useful correspondence. 
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